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Abstract. 

XMM-Newton EPIC observations of the planetary nebula (PN) NGC7009, the Saturn Nebula, have detected 
extended X-ray emission in its central cavity. The diffuse X-ray emission must originate in the shocked fast stellar 
wind. The spectra show that the temperature of the hot gas is 1.8 10 6 K. The rms density derived from the volume 
emission measure is a few tens H-atom cm -3 . The hot gas does not appear over-pressurized with respect to the 
nebular shell. NGC 7009 represents an evolutionary stage at which the influence of the hot gas in the PN interior 
starts to decline due to the diminishing strength of the fast stellar wind and the expansion of the central cavity. 
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1. Introduction 

It has been commonly accepted that the shell structure 
of a planetary nebula (PN) is formed by the current fast 
stellar wind sweeping u p the previo us asymptotic giant 
branch (AGB) wind (e.g. Kwok, 1983 ). In this interacting- 
stellar-winds model, a PN will comprise a central cavity 
filled with shocked fast wind at temperatures of 10 6 — 
10 s K, a dense shell of swept-up AGB wind at 10 4 K, 
and an outer envelope of unperturbed expanding AGB 
wind. The shocked fast wind has a high temperature but 
does not emit appreciably in X-rays because its density 
is low. At the interface between the shocked fast wind 
and the dense nebular shell, heat conduction lowers the 
temperature and mass evaporation raises the density of 
the shocked fast wind, producing optimal conditions for 
X-ray emission ( Weaver et al., 1977 ). Thus, soft X-ray 
emission from the shocked fast wind in a PN interior is 
expected to peak near the inner wall of the nebular shell. 

The ROSAT X-ray satellite had been used to search 
for diffuse X-ray emission from PNe with limited success 
( |Guerrero et al. , 200C ) . It is only recently that the Chandra 
X-ray Observatory has unambiguously detected diffuse X- 
ray emission in the central cavities of three young PNe: 
BP +30° 3639 ( |Kastner et al., 2000|), NGC 65 43 jChu et 
al., 20011) , and NGC 7027 ( [Kastner et al., 200l|) . In the case 
of NGC 6543 where the X-ray emission is well- resolved, a 
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clear limb-brightening is observed. The plasma temper- 
atures implied by their X-ray spectra are a few 10 K. 
Chandra observations of the evolved PN NGC 7293 (the 
Helix Nebula), on the other hand, does not show any dif- 
fuse X-ray emission. These results suggest that the hot gas 
content of a PN diminishes as the PN ages. To confirm this 
hypothesis, PNe at intermediate evolutionary stages need 
to be observed. 

We have examined all pointed and serendipitous ob- 
servations of PNe in the ROSAT archive and found that 
NGC 7009, the Saturn Nebula, hosts an X-ray source with 
spectral characteristics of plasma emission. NGC 7009 is 
more evolved than the three PNe with diffuse X-ray emis- 
sion because of its larger nebular size and dynamical 
age. Consequently, NGC 7009 is an excellent candidate 
not only for diffuse X-ray emission, but also as a PN 
at intermediate evolutionary stages. We have obtained 
X-ray observations of NGC 7009 with the XMM-Newton 
Observatory, and found that the X-ray emission is indeed 
extended. This paper reports our analysis of the diffuse 
X-ray emission from NGC 7009. 

2. Observations 

NGC 7009 was observed with the XMM-Newton 
Observatory in Revolution 255 on 2001 April 30 us- 
ing the EPIC/MOS1, EPIC/MOS2, and EPIC/pn CCD 
cameras. The two EPIC/MOS cameras were operated in 
the Full-Frame Mode for a total exposure time of 39.4 ks, 
while the EPIC /pn camera was operated in the Extended 
Full-Frame Mode for a total exposure time of 31.6 ks. 
The EPIC/pn observations started about 4 ks later 
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than the EPIC/MOS observations because of a larger 
overhead required by the EPIC /pn, but both observations 
concluded at the same time. For all observations, the 
Medium filter was used. At energies below 1 keV, the 
point spread function (PSF) of the EPIC instruments has 
a small core with FWHM ~ 6" but extended wings with 
a half energy width of - 15". For the EPIC/MOS and 
EPIC/pn CCDs, the pixel sizes are l'.'l and 4'.'1, and the 
energy resolutions are ~ 60 eV and ~ 100 eV at 1 keV, 
respectively (Dahlem 1999; Striider et al. 2001). 

We received the XMM pipeline products which in- 
cluded the raw and processed data of all instruments. 
These data were further processed using the XMM- 
Newton Science Analysis Software (SAS version 5.2.0) and 
the calibration files from the Calibration Access Layer 
available on 2002 January 24. 

The event files were further screened to eliminate bad 
events, such as those due to charged particles, and peri- 
ods of high background. For the EPIC/MOS observations, 
only events with CCD patterns 0-12 (similar to ASCA 
grades 0-4) were selected; for the EPIC/pn observation, 
only events with CCD pattern (single pixel events) were 
selected. To assess the background rate, we binned the 
counts over 25 s time intervals for each instrument in the 
10-12 keV energy range, where the counts were dominated 
by the background. The time intervals with high back- 
ground, i.e. count rates > 0.4 cnts s _1 for the EPIC/MOS 
or > 1.5 cnts s _1 for the EPIC/pn, were discarded. The 
resulting exposure times are 35.0 ks, 35.4 ks, and 26.2 ks 
for the EPIC/MOS1, EPIC/MOS2, and EPIC/pn obser- 
vations, respectively. A total of 480±25 cnts, 470±25 cnts, 
and 1, 610 ± 45 cnts from NGC 7009 were detected in the 
EPIC/MOS1, EPIC/MOS2, and EPIC/pn observations, 
respectively. 



3. Spatial Analysis 

We have analyzed the spatial distribution of the X-ray 
emission from NGC 7009 and compared it to the optical 
nebula. As the data analysis is fairly complex, we present 
images at various stages of the processing in Figure 1. 
First, the event files from the three EPIC instruments 
are merged together to construct a raw EPIC image with 
a pixel size of 1"5 in the 0.15-1.25 keV energy band 
(Fig. la). As shown in §4, this energy band includes all 
photons detected from NGC 7009. The raw EPIC image 
is then adaptively smoothed using Gaussian profiles with 
FWHM ranging from 5" to 17". A grey-scale presenta- 
tion of this adaptively smoothed EPIC image is shown in 
Fig. lb. The contours from this smoothed EPIC image 
are overplotted on an Ha image of NGC 7009 taken with 
the Hubble Space Telescope (HST) Wide Field Planetary 
Camera 2 (WFPC2) and presented in Fig. lc. Finally, to 
improve the spatial resolution, this smoothed EPIC image 
is deconvolved using the Lucy- Richardson algorithm and 
a theoretical EPIC PSF provided in the XMM-Newton 
Calibration Access Layer. The contours extracted from 



the PSF-deconvolved EPIC image are overplotted on the 
HST WFPC2 Ha image in Fig. Id. 

The X-ray emission from NGC 7009 is clearly ex- 
tended. Its spatial distribution in the adaptively smoothed 
EPIC image (Fig. lb and lc) can be described by an el- 
liptical Gaussian with FWHM ~ 23" x 14" and a major 
axis along the position angle ~ 70°. There is an excel- 
lent correspondence between the X-ray emission and the 
nebular shell of NGC 7009. The ffi>TWFPC2 Ha image of 
NGC 7009 shows a bright 25" x 10" inner shell surrounded 
by a 25" x 20" envelope, with two ansae extending from 
the tips of the major axes of the inner elliptical shell out 
to 25" from the central star. The position and orientation 
of the 50% contour of the smoothed EPIC image matches 
almost perfectly the inner shell of the nebula (Fig. lc), in- 
dicating that the X-ray emission peaks within the central 
cavity. The PSF-deconvolved EPIC image further shows 
that the majority of the X-ray emission originates from the 
central cavity. As the EPIC PSF has extended wings and 
our PSF-deconvolution is limited by the number of counts 
and the validity of the theoretical PSF adopted, it is most 
likely that all diffuse X-ray emission from NGC 7009 is 
confined within the inner shell, similar to what is observed 
in BD+30°3639 and NGC 6543. Furthermore, NGC 7009 
has a sharp-edged inner shell similar to that of NGC 6543, 
which implies an expansion driven by the pressure of the 
hot gas in the central cavity (Chu et al., 2001). Therefore, 
we suggest that the hot, X-ray-emitting gas in NGC 7009 
is also confined within its central cavity, but the detailed 
distribution of the hot gas cannot be determined from the 
currently available X-ray observations. In particular, it is 
not possible to address whether X-ray emission from its 
central star is present or not. 

4. Spectral Analysis 

Spectra of NGC 7009 have been extracted from the event 
files of the three XMM-Newton EPIC cameras separately 
because they have different instrumental response func- 
tions. A 50" x 40" elliptical source aperture is used to 
include all X-ray emission detected from NGC 7009. An el- 
liptical annular region exterior to and ~ 10 times as large 
as the source aperture is used to assess the background 
level. The background-subtracted spectra of NGC 7009 are 
shown in Fig. 2- left. 

The EPIC/pn spectrum peaks at ~ 0.55 keV and 
drops abruptly to a faint plateau between 0.7 and 1.0 keV, 
but stays essentially flat from 0.5 keV down to 0.3 keV. 
The EPIC/MOS spectra also peak at ~ 0.55 keV, but 
have smaller numbers of counts because of lower sensitiv- 
ities. The detailed differences between the EPIC/pn and 
EPIC/MOS spectra reflect the different instrumental re- 
sponses among these cameras and the low S/N ratios in 
the latter spectra. Neither the EPIC/MOS spectra nor 
the EPIC/pn spectrum shows any detectable emission at 
energies greater than 1.25 keV. 

The EPIC spectra of NGC 7009 are clearly inconsis- 
tent with a power-law or blackbody spectrum. Thus, the 
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X-ray emission from NGC 7009 must be thermal plasma 
emission and the spectral peak at ~ 0.55 keV corre- 
sponds to the He-like triplet of O vn at ~ 0.57 keV. 



gas, 1.8 10 6 K, is below that expected for the shocked fast 
wind from its central star which has a terminal velocity 
~ 2,800 km s" 1 (|Cerruti-Sola fc Pcrinotto, 1989|; [fping 



We car. model the observed spectra in order to determine & Sonneborn, 2001). Similar situation is observed in the 



the physical conditions of the X-ray-emitting plasma and 
the amount of foreground absorption. We have adopted 
the thin plasma emission model of Raymond fc Smitlj 
(1977|yand the nebular abundances determined from opti- 



other PNe with diffuse X-ray emission: the plasma tem- 
peratures implied by the X-ray spectra of BD +30°3639 
and NGC 6543 are several times lower than the expected 



post-shock temperature for their fast stellar wind (Kastnci 



cal anc UV spectra of NGC 7009. The abundances of He, 
C, N, O, Ne, and S relative to hydrogen by number have 
been measured to be 0.13, 2.9 10~ 4 , 2.2 10~ 4 , 4.3 10" 4 , 
7.2 10~ 5 , and 1.9 10 -5 , respectively (Kwitter & Henry 
1998). Accordingly, we have a dopted nebular abundances , 
relative to the solar values ( Anders fc Grevesse, 1989] ), 
of 1.33, 0.8, 1.9, 0.7, 0.6, and 1.2 for He, C, N, O, Ne, 
and S, respectively, and 1.0 for the other elements. We 
have assumed solar abundances for the foreground inter- 
stellar absorption, and adopted absorption cross-sections 
from Balucinska-Church fc McCammon (1992 ). 



et al., 200C; Chu et al., 2001). This temperature difference 



The spectral fits are carried out by folding the 
Raymond-Smith model spectrum through the correspond- 
ing EPIC/MOS1, EPIC/MOS2, and EPIC/pn response 
matrices, and fitting the modeled spectra to the observed 
EPIC spectra in the 0.3-1.5 keV energy range using the % 2 
statistics. The best- fit model, overplotted on the EPIC/pn 
and EPIC/MOS spectra in Fig. 2-left, has a plasma tem- 
perature of T = 1.8 10 6 K (or kT = 0.152 keV), an ab- 
sorption column density of TVh = 8.1 10 20 cm~ 2 , and a 
volume emission measure of 1.9 10 54 gP cm~ 3 where d is 
the distance to NGC 7009 in kpcQ. 

The goodness of the spectral fits is illustrated by the 
plot of reduced \ 2 of the fits as a function of kT and 
A^h shown in Fig. 2-right. The 99% confidence contour 
spans T = 1.6-2.0 10 6 K (or kT = 0.14-0.17 keV), A H 
= 2—17 1 20 cm -2 , and volume emission measure 1.5- 
2.4 1 54 d 2 cm -3 . Assuming that the HI column density is 
similar to the absorption column density and adopting a 
typical gas-to-dust ratio (Bohlin, Savage, & Drake 1978), 
the co rresponding visual extinction is Ay = 0.44^° ^ mag, 



is expected if heat conduction is prevalent at the interface 
between the shocked fast wind and the dense nebular shell 
( |Soker, 1994 |Zhckov fc Pcrinotto, 1996] ). 

In order to determine the role played by the hot 
gas in the nebular dynamics of NGC 7009, it is possi- 
ble to derive the rms density and pressure of the X- 
ray-emitting gas and compare them to those of the sur- 
rounding cool nebular shell. Assuming a prolate ellip- 
soidal central cavity, the volume occupied by the X-ray- 
emitting gas is ~6.7 10 51 ((i/1.2) 3 e cm 3 , where e is the 
filling factor. From this volume and the aforementioned 
volume emission measure we derive an rms density of 
20 (d/1.2)- 1 / 2 e^ 1 / 2 cm" 3 , and a total mass for the X- 
ray-emitting gas of 8 10~ 5 (d/1.2) 5 / 2 e 1 / 2 M Q . This rms 
density and the plasma temperature of 1.8 10 6 K imply 
that the thermal pressure of the X-ray-emitting gas is 
-4 10~ 9 (d/1.2)- 1 / 2 e" 1 / 2 dynes cm" 2 . The nebular shell 
surrounding the hot gas has a density of 5,000 cm~ 3 , a 
temperature of 10,500 K (Hyung & Aller 1995; Rubin et al. 
2001), and thus a thermal pressure of ~T0~ 8 dynes cm" 2 . 
The thermal pressure of the nebular shell is comparable 
or higher than that of the X-ray-emitting gas unless the 
filling factor e is <C0.2. 

The X-ray luminosity of NGC 7009 and the density and 
pressure of the hot gas in its interior are the lowest among 
the four PNe for wich diffuse X-ray emission has been de- 
tected thus far. These differences may be associated with 
the nebular evolution. NGC 7009 has a large nebular shell, 
0.145x0.06 pc in size, and dynamical age, ~ 1, 700 yrs for 
an expansion velocity of 17 km s" 1 ( Bohigas, Lopez, fc 



and the logarithmic extinction at the H/3 line is egg — Aguilar, 1994). Therefore, it appears to be in the most 



0.15 



+0^7 
-0.15- 



This extinction is consistent with other esti- 



mates of extinction (ch/3 = 0.24 ± 0.06, Lame & Pogge 
1996; ch/3 = 0.10+o EJg, Rubin et al. 2001). The observed 
(absorbed) X-ray flux of NGC 7009 in the 0.3-1.5 keV en- 
ergy range is 7.2 10~ 14 ergs cm -2 s _1 , and the unabsorbed 
(intrinsic) X-ray flux is 1.7 10~ 13 ergs cm~ 2 s _1 . The total 
X-ray luminosity in the 0.3-1.5 keV energy range is Lx = 
2.0 10 31 d 2 ergs s~ 4 = 3 10 31 ergs s" 1 for a distance of 1.2 
kpc. 



5. Discussion 

XMM-Newton EPIC observations of NGC 7009 have 
shown that the diffuse X-ray emission is confined within 
the central cavity. The temperature of the X-ray-emitting 



The distance estimates to NGC 7009 range from 0.5 to 2.1 
kpc. In the following we will adopt a distance of 1.2 kpc. 



advanced evolutionary stage among these four PNe. At 
a more advanced evolutionary stage, NGC 7293 does not 
have any detectable diffuse X-ray emission (Guerrero et 
al. 2001). The Chandra and XMM-Newton observations 
of these five PNe suggest that PNe are brightest in X-rays 
when they are young and NGC 7009 may be at an evo- 
lutionary stage when the diffuse X-ray emission starts to 
decline. This evolutionary effect is expected as the stel- 
lar wind diminishes and the central cavity grows. More 
observations of PNe at different evolutionary stages are 
needed to evaluate the duration of the pressure-driven 
phase in PNe and its effects on the nebular dynamics. 
As the evolved PNe are expected to have low X-ray sur- 
face brightness, the high sensitivity of XMM-Newton is 
essential in detecting such PNe. 
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Fig. 1. Raw (a) and smoothed (b) XMM- Newton EPIC images of NGC 7009 in the 0.15-1.25 keV energy band. X-ray 
contours from the smoothed EPIC image (c) and from the PSF-deconvolved EPIC image (d) overplotted on an HST 
WFPC2 Ha image of NGC 7009. The contour levels are chosen at 10, 25, 50, 75, and 90% of the emission peak. 




Energy [keV] kT [keV] 

Fig. 2. (left) XMM-Newton EPIC/pn (top), and EPIC/MOS1 and EPIC/MOS2 (bottom) background-substracted 
spectra of NGC 7009. The spectra are overplotted with the best-fit Raymond & Smith model (1977) to the EPIC/pn 
and the combined EPIC/MOS spectra. Both spectra and models are plotted in 30 eV wide energy bins to preserve 
the EPIC spectral resolution (~ 60-100 eV). (right) x 2 grid plot of the spectral fit of the EPIC/pn and the combined 
EPIC/MOS spectra. The contours represent 68%, 90%, and 99% confidence levels, and the '+' sign corresponds to 
the parameters of the best-fit model. 



